Thermodynamic properties of a magnetically modulated graphene 
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The effect of magnetic modulation on thermodynamic properties of a graphene monolayer in pres- 
ence of a constant perpendicular magnetic field is reported here. One-dimensional spatial electric 
or magnetic modulation lifts the degeneracy of the Landau levels and converts into bands and their 
bandwidth oscillates with magnetic field leading to Weiss-type oscillation in the thermodynamic 
properties. The effect of magnetic modulation on thermodynamic properties of a graphene sheet is 
studied and then compared with electrically modulated graphene and magnetically modulated con- 
ventional two-dimensional electron gas (2DEG). We observe Weiss-type and de Haas- van Alphen 
(dHvA) oscillations at low and high magnetic field, respectively. There is a definite phase differ- 
ence in Weiss-type oscillations in thermodynamic quantities of magnetically modulated graphene 
compared to electrically modulated graphene. On the other hand, the phase remains same and am- 
plitude of the oscillation is large when compared with the magnetically modulated 2DEG. Explicit 
asymptotic expressions of density of states and the Helmholtz free energy are provided to understand 
the phase and amplitude of the Weiss-type oscillations qualitatively. We also study thermodynamic 
properties when both electric and magnetic modulations are present. The Weiss-type oscillations 
still exist when the modulations are out-of-phase. 

PACS numbers: 65.80.+n,71.70.Di,71.18.+y 



I. INTRODUCTION 

Graphene is a two-dimensional sheet of carbon al- 
lotrope with honeycomb lattice structure. It can be con- 
sidered as a basic building block of all other dimensional 
carbon allotropei~— . The relativistic-like, massless and 
linear energy spectrum of graphene's quasi-particles in 
low-energy range close to the Dirac points in it's band 
structure reflect in it's different properties like transport 
properties, optical properties, magnetic properties etc in 
different way than the conventional 2DEG formed in the 
semiconductor heterostructures. The massless linear en- 
ergy dispersion and charge conjugation symmetry cause 
some unusual phenomena like Klein paradox, anomalous 
quantum Hall effect and non-zero Berry phased—. 

Effect of electric or/and magnetic modulations on 
transport and thermodynamic properties of quantum 
2DEG systems is continuing to be an active research 
field. The magnetotransport properties of a conventional 
2DEG in presence of a one-dimensional (ID) weak elec- 
tric modulation have been studied in great details ex- 
perimentally and theoretically by various groups^—. 
Later, the magnetotransport properties of a magnetically 
modulated 2DEG in presence of a constant perpendic- 
ular magnetic field have been studied theoretically-i 6 ^— 
and also experimentally-- - — . The presence of weak elec- 
tric/magnetic modulation broadens the Landau energy 
levels into bands. The band width oscillates with the 
magnetic field and its oscillations are also reflected in 
magnetotransport properties. It has been observed that 
the magnetoresistivity tensor oscillates with inverse of 
the magnetic field at very low magnetic field. This os- 
cillation is known as the Weiss oscillation which is com- 
pletely different from the Shubnikov-de Hass (SdH) os- 
cillations observed at higher magnetic field. Period of 



Weiss oscillation varies with electron density n e as i/n e , 
whereas for SdH it depends linearly on n e . The Weiss os- 
cillation is due to effect of the commensurability between 
two length scales in the system: the cyclotron diameter 
at the Fermi energy and the period of the modulated elec- 
tric/magnetic potential. Magnetotransport properties of 
electricall y 23 ' 24 and magnetically modulated graphene 2 ^ 
in presence of a constant perpendicular magnetic field 
have been studied recently. 

The effect of weak electric and magnetic modula- 
tion on thermodynamic properties for 2DEG in pres- 
ence of a perpendicular magnetic field have been stud- 
ied theoreticall y 12 ' 26 ' 27 . It is observed that the Weiss- 
type oscillation in various thermodynamic properties are 
present. The Weiss- type oscillations are completely dif- 
ferent in origin from the usual de Haas-van Alphen 
(dHvA) oscillations which occurs at high magnetic field. 
The dHvA oscillation is effect of the formation of dis- 
crete Landau energy levels due to the quantizing mag- 
netic field. Recently, thermodynamic properties of a 
monolayer graphene in presence of weak electric modu- 
lation have been studied and the Weiss-type oscillations 
are seen 2 ^. These results inspired us to study thermody- 
namic properties of a magnetically modulated graphene 
sheet in presence of a constant magnetic field. 

The source of the magnetic modulation, magnetic 
stripes or superconducting materials, acts like electrical 
gates and produces an electric modulation. The trans- 
port properties of a 2DEG 16 and a graphene sheets in 
presence of both the modulations were studied. 

In this paper we study the effect of magnetic modula- 
tion on thermodynamic properties of a graphene sheet in 
presence of a constant magnetic field. We compare our re- 
sults with electrically modulated graphene and with mag- 
netically modulated conventional 2DEG. We also calcu- 
late an asymptotic expression of density of states (DOS) 



and the Helmholtz free energy of a magnetically modu- 
lated graphene in presence of a constant magnetic field 
at low temperature. In addition to these, we also study 
thermodynamic properties of graphene and 2DEG when 
both electric and magnetic modulations are present. 

This paper is organized as follows. In section II, 
we summarize the standard results of the energy spec- 
trums and the corresponding eigenstates for electrically 
and magnetically modulated graphene layer and 2DEG 
in presence of a constant magnetic field. In section 

III, we numerically calculate the thermodynamic quanti- 
ties like, Helmholtz free energy, internal energy, entropy, 
heat capacity and magnetization. We analyze our nu- 
merical results and compare with electrically modulated 
graphene and magnetically modulated 2DEG in section 

IV. We discuss the behavior of the asymptotic expression 
of the DOS and the Helmholtz free energy for magneti- 
cally modulated graphene in section V. In section VI, we 
study thermodynamic properties of graphene and 2DEG 
in presence of both electric and magnetic modulations. 
We present summary of our work in section VII. The 
detail calculation of the DOS by using Green's function 
method is presented in the Appendix 1. 



II. ENERGY SPECTRUM 

We consider a monolayer graphene sheet subjected to 
a perpendicular constant magnetic field Bo = Boz, the 
Hamiltonian of an electron with charge — e is given by 



*? 



v F a ■ (p + eA ) 



(1) 
10 6 



where er = (a x ,a y ) are the Pauli matrices, v F 
m/s is the Fermi velocity and Ao is the magnetic vec- 
tor potential. Here, we have chosen the Landau gauge 



Ao = Boxy. The energy eigenvalues are E% = hui g \j2n 1 
where n = 0,1, 2,3... is the Landau level index and 
u! g = v F y/eBo/h is the cyclotron frequency. The cor- 
responding normalized eigenstates are 



^n,k y (x,y) 



ik y y 



y/2Lylo 



-i<t>n-l[(x + Xo)/lo] 
<f>n[(x + Xq)/1 ] 



(2) 



where 4>n{x) is the known harmonic oscillator wave func- 
tion, l = y/h/(eBo) is the magnetic length, x = — k y l 2 
is the center of the cyclotron orbit and L y is the width 
of the graphene in the y-direction. 

We consider the perpendicular magnetic field is modu- 
lated very weakly by Bi = B\ cos (qx)z, where q = 2n/a 
and a is the modulation period. The total Hamilto- 
nian can be written as H = H$ + Hi, where Hi — 
V^a y sin (qx). Here, with V^ = ev F Bi/q is the strength 
of the effective magnetic potential determined by the am- 
plitude Bi and the period a of the magnetic modulation. 
Using the first-order perturbation theory, the energy cor- 
rection to the unperturbed Hamiltonian Hq is given as2^ 

A£ S" = V^e^ 2 [L n . 1 (u)-L n (u)]cos{qx %3) 



where L n (u) is the Laguerre polynomial and u — q 2 l\j2. 
So the total energy upto the first-order in V^ is given by 
E%™ = El + AE g n ™ y . The band width in presence of the 

magnetic modulation is A TO ~ sin(2y / ml — 7r/4). Using 
the flat band condition, we get 2R C = a(j + 1/4), with 
j = 0,1,2,3,... andi? c = k F l 2 . 

The energy correction due to the weak electric modu- 
lation U(x) = U/ cos (qx) can be obtained by the same 
method and it is given as 2 ^ 

V 9 _ /9 

-%-e u ' z [L n (u) + L n -i(«)] cos(qx ). 






(4) 



So the total energy upto the order of V£ is E^' e k — 
E% + AE 9 't . The bandwidth in presence of the electric 
modulation is A e ~ cos(2y / rm — 7r/4). The bandwidths 
due to magnetic modulation and due to the electric mod- 
ulation are out of phase. The band will be flat when 
2R C = a(j + 3/4) which is different from what we get in 
the magnetic modulation case. 

The Hamiltonian of a conventional 2DEG in presence 
of a perpendicular constant magnetic field Bq is 
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The energy spectrum is E^ d = huj c (n + 1/2), where n = 
0, 1, 2, 3.... and u; = eB /m* is the cyclotron frequency. 
The corresponding eigenstates are 
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In presence of the weak magnetic modulation Bi, the 
total Hamiltonian is if = H]f + Hf d , where H\ d — 
[V£ d /(hq)](p y + eB x)sm(qx) with V™ = H(eBi/m*) 
is the strength of the effective magnetic potential de- 
termined by the amplitude of the magnetic modulation. 
The first-order energy correction due to the weak mag- 
netic modulation is given b y 16 i 



AE; 



2ri.i 
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The total energy of a 2DEG in presence of the magnetic 
modulation is then E 2d ^ = £ 2rf + AE 2 n d j^ . The energy 
correction due to electric modulation is given by 

y2d 
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u)]cos(qx ). (8) 



All these standard results will be used to calculate the 
thermodynamic properties numerically in the next sec- 
tion. 



III. THERMODYNAMIC QUANTITIES 

In this section we discuss all standard thermodynamic 
equations to be used for calculating chemical potential, 



Helmholtz free energy, internal energy, entropy, magne- 
tization and specific heat. 

The DOS of a magnetically modulated graphene sheet 
in presence of a constant magnetic field can be written 
as 



D(E) 



A 



E^ 



E„ 



(9) 



where A = L x L y is the area of the graphene sheet 



and E Ui k is a energy dispersion of a given system like 
^nk i Ent ' E n j™ 1 . The dependence of chemical poten- 
tial fi(B,T) on temperature (T) and magnetic field (B) 
can be obtained numerically by using the following nor- 
malization condition 



N 



D(E)f(E)dE, 



(10) 



where N is the total number of electrons, f(E) — 
[exp(jr-r) + 1] _1 is the Fermi-Dirac distribution function 
and k B is the Boltzmann constant. Using the expression 
of D(E) given in the above equation we get 



^0 = " E / dt f( E n,t), 



(11) 



where n e is the electron density and t = qxQ. 
The total internal energy can be written as 



U= / ED(E)f(E)dE. 
Jo 

The internal energy per unit area is 



U I 
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- 271 /, / E n> tf (E„ it )dt. 
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(12) 



(13) 



Now for a system of non-interacting electrons, the 
Helmholtz free energy density 30 is given by 



k a T 



-r = me 
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n=0 
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En 



k a T 



In 1 + exp ( — ; — ^- L - ) dt. 

(14) 

From the above equations it is clear that the DOS plays 
an important role in the behavior of the thermodynamic 
properties. In presence of the perpendicular magnetic 
field the DOS shows a series of delta function because of 
the quantized energy spectrum. Graphene and conven- 
tional 2DEG having different energy spectrums reflects 
differently in their thermodynamic properties. By using 
the above results we compute entropy via S = (U—F)/T, 
orbital magnetization M = — (dF/dBo)A,N and heat ca- 
pacity C = T{dS/dT) A ^ N = -T{d 2 F/d 2 B ) A ,N- For 
better visualization of effect of the magnetic modulation 
we plot the fluctuation All = II(Bi) - U(B X = 0), where 
Bi is the strength of magnetic modulation, and IT is a 
thermodynamic quantity like /i, F, U, M, S, C . 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

Thermodynamic properties of magnetically modulated 
graphene sheet in presence of a constant magnetic field 
are studied. The aim is to study the effect of magnetic 
modulation on graphene in comparison with the electri- 
cally modulated graphene and magnetically modulated 
conventional 2DEG. We plot the fluctuation due to weak 
modulation in chemical potential, Helmholtz free energy, 
magnetization, internal energy, entropy and specific heat 
with the magnetic field B n . We have used the follow- 
ing parameters for numerical calculation: electron den- 
sity n e — 3.16 xlO 15 / m 2 , effective mass of an electron 
rn* = 0.067m e with m e is the bare electron mass, tem- 
perature T = 2 K, modulation period a — 382 nm as 
used i n 26 i 28 . For these parameters, E 9 F = 0.1 eV and 
Ef = 14.25 meV. 




FIG. 1: Plots of the fluctuation in chemical potential vs mag- 
netic field at T = 2 K. 

In figures 1-6, we have plotted the fluctuations in var- 
ious thermodynamic properties, An, due to both mag- 
netic and electric modulations. We have scaled all the 
thermodynamic quantities per electron. Each figure con- 
tains two panels, the upper panel shows the effect of the 
magnetic modulation on graphene (solid line) and con- 
ventional 2DEG (dashed line) , and the lower panel shows 
this fluctuation for the electrically modulated (dashed 
line) and magnetically modulated (solid line) graphene. 
We have zoomed the oscillations at low magnetic field 
and shown in the inset of figures 1-6. It clearly shows 
that weak ID periodic potential, either electric or mag- 
netic in nature, induces new oscillations at low magnetic 




FIG. 2: Plots of the change in the free energy vs magnetic 
field. To make the fluctuation dimensionless we use Fq — 
Ep/2 in the lower panel but as Ef is different for the two 
systems we are showing fluctuation in units of k B T in the 
upper panel. 



field. These modulation induced oscillation is due to the 
commensurability of the two length scales present in the 
system. These oscillations are similar to the Weiss os- 
cillations observed in the magnetoresistance at low mag- 
netic fields. 

It is clear from the upper panels of the figures 1-6 
that the fluctuation in the thermodynamic properties, 
All, of the graphene sheet is quite large compared to 
the 2DEG system. This can be understood qualitatively 
from the following arguments. The energy correction for 
the magnetically modulated graphene and 2DEG sys- 
tems are VJ^ = 1 meV and V%f = 0.046 meV for the 
same strength of the magnetic modulation B\ — 0.02 T. 
Clearly, the energy correction due to the magnetic mod- 
ulation in graphene is quite large compared to that of the 
2DEG. This is the origin for higher amplitude fluctuation 
of the Weiss-type and dHvA oscillations in graphene com- 
pared to the 2DEG. It is interesting to note that ampli- 
tude of diffusive conductivity in magnetically modulated 
graphene is small compared to magnetically modulated 
2DEG2£. 

The lower panels of the figures 1-6 show that the Weiss- 
type oscillation has a definite phase difference in the 
thermodynamic properties between magnetic and elec- 
tric modulation cases, which is due to the following rea- 
sons. The phase difference in the fluctuations between 
electrically and magnetically modulated graphene comes 



from the nature of the energy correction. In the electric 
modulation case the energy correction contains addition 
of the two successive Laguerre polynomials and in mag- 
netic modulation case it is subtraction of the two suc- 
cessive Laguerre polynomials, giving rise to the cosine 
and sine term, respectively. We also observe that the 
dHvA-type oscillations for the two different kind of mod- 
ulations remain in the same phase with each other. The 
phase difference between the two different kind of modu- 
lated systems are shown explicitly by using the analytical 
expressions of the DOS in the next section. 

Figures (1) and (3) shows the fluctuation in chemi- 
cal potential and magnetization with magnetic field, re- 
spectively. The fluctuation in magnetically modulated 
graphene is several times higher than the magnetically 
modulated 2DEG in the Weiss-type oscillation, but their 
phases remain same. On the other hand, the magnet- 
ically modulated graphene shows 7r phase difference in 
Weiss-type oscillation compared to the electrically mod- 
ulated graphene with the same amplitude. 

Figure (2) is showing the fluctuation in the Hclmholtz 
free energy where lower panel shows tt/2 phase differ- 
ences between electric and magnetic modulation cases 
in graphene with the same amplitude. The amplitude 
of fluctuation of graphene is several times higher than 
the 2DEG but the phases remain same. From the flat- 
band condition, the minima of the bandwidth occur at 
B(T) = 0.092,0.113,0.148,0.214,0.386. On the other 
hand, the fluctuation in the free energy vanishes at 
B{T) = 0.092,0.115,0.150,0.214,0.385. It shows that 
the minima of the free energy fluctuation occur at those 
values of the magnetic field where the bandwidth minima 
occur. 

We have plotted internal energy fluctuation in Fig. (4) . 
The lower panel of Fig. (4) shows that the Weiss-type os- 
cillation is appeared with same amplitude but tt/2 phase 
difference when compared with the electrical modulation 
case. When we compare with conventional 2DEG, it is 
similar to the case of the Hclmholtz free energy i. e. 
amplitude is higher in magnitude. 

Figures (5) and (6) are showing the entropy and spe- 
cific heat fluctuation, respectively. The fluctuation in 
graphene is higher by several times than that of the 
2DEG. The phase relationship is not discernible in these 
figures. 

In all the above cases, dHvA-type oscillation remains 
in the same phase and does not depend on the modula- 
tion type as it is the manifestation of the quantized Lan- 
dau levels rather than periodic perturbation. Figure (7) 
shows the damping of fluctuation in chemical potential 
and free energy with increasing temperature. Tempera- 
ture dependence in Weiss-type and dHvA oscillations are 
independent of type of modulation and already discussed 
in electrical modulation case. The phase relationships for 
all the thermodynamic quantities are given in Table (1). 

Even at higher magnetic field, the effect of modulation 
on the fluctuation of thermodynamic quantities is still 
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FIG. 3: Plots of the fluctuation in magnetization vs magnetic 
field. In the lower panel AM is scaled by Mo = fi* B , where 
fx B — eh/2m g is the effective Bohr magneton with m g — 
E 3 F /v 2 F is the cyclotron mass. 



Case 


A^i 


AF 


AM 


AU 


AS 


AC 


Phase shift 


7T 


tt/2 


7T 


tt/2 


indiscernible 


indiscernible 



TABLE I: The phase shifts in the Weiss-type oscillations ap- 
pear in the fluctuation of thermodynamic quantities between 
the electrically and magnetically modulated graphene. 



exist which is the manifestation of the modulated density 
of states. The dHvA-type oscillation is corresponding to 
the crossing of each Landau level one by one through the 
Fermi level. The density of available states per Landau 
level is 1/2ttIq and it increases linearly with magnetic field 
which results in increasing amplitude of dHvA oscillation 
with magnetic field. 



V. ASYMPTOTIC RESULTS 



FIG. 4: Plots of the fluctuation in the internal energy vs 
magnetic field. In the lower panel AU is scaled by Uo = Ep/2. 



An(e) = - 



2fL 



7TZq HUJg 



(7T 
^° e — 7 



and 



(16) 



,2 , „\ 3 / l \ ■• 



^=^-(t) (^-) -Mt-t)- iv.-) 



IT \Iq J \ tkd 



2e 9 



Here, e = E/{huj g ) and e% = E 9 F /(hu g ). 

The fluctuation in the DOS for electrically modulated 
graphene^ is proportional to cos 2 (qloe 9 F — 7r/4). The ap- 
pearance of the square of the sine term in Eq. (| 16|) in- 
stead of square of cosine is the reason of definite phase 
differences in the fluctuation of all the thermodynamic 
quantities. 

Using the two Eqs. (fl"5*)) and (fT()|) separately, we get 
an approximate analytical expression of the free energy 
and it's fluctuation. Our aim is to study the magnetic 
modulation effect in compare to the electrical modulation 
in graphene^. The change in the free energy due to the 
magnetic modulation can be expressed as 



Here we derive an asymptotic expression for the DOS 
and the Helmholtz free energy. For weak magnetic mod- 
ulation and under quasi-classical limit, we calculate the 
DOS by using the Green's function technique (see the 
Appendix 1) and written as a sum of the modulated and 
unmodulated part as D(e) = D u (e) + Z) m (e), where 
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TTIq fVjJg 



1 + 2 cos (?re 2 
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1 + exp 



k a T 



de. 



dc 



Under the assumption of very low temperature the above 





FIG. 5: Plots of the fluctuation in entropy vs magnetic field. 
In the lower panel AS is scaled by So = k B . 



FIG. 6: Plots of the fluctuation in specific heat vs magnetic 
field. 
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where T g iHvA 



(faujg)/(2ir 2 k B e 9 F ) is the critical temper- 
ature for the dHvA-type oscillations in graphene. The 
ratio of amplitude of the free energy fluctuations of the 
magnetically and electrically modulated graphene is 



AI 



\vi) ■ 



(19) 



Here, A m and Af are the amplitudes of free energy 
fluctuation for magnetically and electrically modulated 
graphene, respectively. The expression of Af is taken 
from the Ref<2£i. But in the case of conventional 2DEG 
this ratio has been calculated in Ref. 26 and it is given by 



2d 



2.1 
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2^2 



Ef 



U 2rf 



(20) 



where e a = h 2 /(m*a 2 ), 7^ and 7 2d are amplitudes of the 
free energy fluctuation for magnetic and electric modula- 
tion, respectively. The ratio of the amplitudes of differ- 
ent modulation cases in graphene behave differently from 
that of the 2DEG system. 



Now we compare our result with the conventional 
2DEG modulated magnetically. In this two case the ratio 
of amplitudes of the free energy fluctuations is given by 



(21) 



where E 9 F is the Fermi energy of graphene whereas Ep d 
is for conventional 2DEG. Using B = 0.1 T, T = 2 K, 
a = 382 nm, f3E F = 521, ujg/ujQ = 35, we get 
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(22) 



This equation leads to X 9 m - 4.7A 2d when B x = 0.02 T. 
In graphene, the amplitude of free energy fluctuation in 
the Weiss-type oscillation is higher and in the same phase 
in compare to the conventional 2DEG. 



VI. ELECTRIC AND MAGNETIC 
MODULATIONS 

In this section, we study how the thermodynamic prop- 
erties discussed in the previous sections are changed in 
presence of an additional electric modulation with the 
same period. We consider two different cases: when 
both the modulations are in-phase and that are 7r/2 out- 
of-phase. The fluctuation in thermodynamic quantities 
like chemical potential and the Helmholtz free energy are 
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FIG. 7: Plots of the fluctuation in chemical potential and free 
energy vs magnetic field for T = 2 K and T = 6 K with the 
same parameter. 



calculated numerically. Other thermodynamic quantities 
can easily be obtained by taking a suitable numerical 
derivative of the Helmholtz free energy fluctuation. 

In-phase modulations: We consider a weak electric 
modulation described by the periodic potential U(x) = 
V£ cos(qx), which is in-phase with the magnetic modula- 
tion Bi(ar) = B\ cos(qx)z. When both the modulations 
are in the same phase, the total energy correction for 
graphenc in a weak magnetic field can be written as 



Ai? 9 - v7£+/!sin(2V^J - - + S?) cos(qx ), (23) 



where f m = (2V9 l / y ^)(n/u 3 ) l / i sm(^/^), f e = 

(KV^)(l/^) 1/4 cos(V^) and Sf = tan" 1 ^//^. 
The flat-band condition at the Fermi energy gives the 
positions of the minima in the free energy fluctuation as 
Bj = 2 P F B a/(j + j--^)- Here, j is an integer, p F = ak F 
is a dimensionless momentum and B a = h/(ea 2 ) is the 
characteristic magnetic field. In this case, Sf — tt/4 and 
then Bj = 2p F B a /j. 

Similarly for conventional 2DEG, the total energy cor- 
rection in the low magnetic field can be written as 



AE 2d ~ y""™ + » e 2 sin(2V^ - - + Sf d ) cos(qx ), (24) 



where w n 
and 5 2d 



tan 1 (w m /w e ). The flat-band condition 



at the Fermi energy is now Bj = 2p F B a /(j + 1/4 — 

tan- 1 {V*'p F /(27rK M )}). 

In figures 8 and 9, we plot the chemical potential 
and free energy fluctuations in presence of both elec- 
tric and magnetic modulations, respectively. The upper 
panel shows the thermodynamic fluctuations of graphene 
(solid) and 2DEG (dashed) when the modulations are in- 
phase and the lower panel shows the thermodynamic fluc- 
tuations of graphene (solid) and 2DEG (dashed) when 
the modulations are out-of-phase. 




FIG. 8: Plots of the fluctuation in chemical potential vs mag- 
netic field. 



Out-of-phase modulation: We consider the same 
electric modulation U(x) = V e cos(qx) and assume mag- 
netic modulation is given by Bi(x) = B\ sin(qx)z so that 
the two modulations are ir/2 out-of-phase. To first-order 
in V e and B±, the total energy correction for graphene 
can be written as 



A£ s = / m Jl + {(^) 2 -l}cos 2 (2V^-l) 

y Jm 

X sm(qx + 5%), 



(25) 



where tan(<5f) = (f m /fe) tan(2 v / ml— 7r/4). Similarly, the 
total energy correction for 2DEG is written as 



AE 



2d 



1 + {(^-)2 - l}cos 2 (2V^ - 1) 



x sin(ga;o + S 2d ) 



(26) 



where ta,n(S 2d ) = (w m /w e ) tan(2 v / mZ — tt/4). In the 
case of graphene, though the band width is almost 
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constant at the Fermi energy, but the magnetic field 
dependent phase factor 5% plays an important role in 
the fluctuation of the thermodynamic quantities. This 
phase factor produces the Weiss-type oscillations in the 
thermodynamic properties at low magnetic field. In 
conventional 2DEG, the band width oscillates with the 
magnetic field when w e =^ w m in our case but the Weiss- 
type oscillation is due to the magnetic field dependent 
phase factor 5% in the total energy correction. When 
w e = w m , the band width becomes non-oscillatory, but 
the Weiss-type oscillation still exist due to the magnetic 
field dependent phase factor. It is interesting to contrast 
our result with the results of the Weiss oscillations in the 
conductivity in presence of both the modulation o 16 ' 29 . 
The Weiss oscillation in conductivity is suppressed when 
the modulations are out-of-phase. On the other hand, 
the Weiss-type oscillations in the thermodynamic prop- 
erties enhanced when the modulations are out-of-phase. 
In out-of-phase case, the Weiss-type oscillation in the 
thermodynamic quantities is due to the magnetic field 
dependent phase factor in the energy correction. 




suits of magnetically modulated graphene are compared 
with electrically modulated graphene and magnetically 
modulated conventional 2DEG. It is observed that in 
the case of magnetically modulated graphene, a definite 
phase difference is appeared in the Weiss-type oscillation 
in compare to the electrically modulated graphene for all 
thermodynamic quantities. But when we compare our re- 
sults with magnetically modulated conventional 2DEG, 
the amplitude of the fluctuations is found to be higher in 
graphene than 2DEG, but the phases remain same. We 
calculate the DOS and the Hclmholtz free energy analyt- 
ically, and explain the origin of the this phase difference 
in the Weiss-type oscillation. The enhancement of the 
fluctuation in magnetically modulated graphene in com- 
pare to the 2DEG and the definite phase difference be- 
tween magnetically and electrically modulated graphene 
are explained by using the asymptotic results of the DOS 
in general and the Hclmholtz free energy in particular. 

We have also studied the thermodynamic properties 
like chemical potential and the Hclmholtz free energy 
of graphene and conventional 2DEG in presence of both 
magnetic and electric modulations with the same period. 
The combined effect of both modulation does not modify 
the nature of the Weiss-type oscillation when they are 
in the same phase except a finite phase shift in the fluc- 
tuation. The effect of the out-of-phase modulations on 
thermodynamic fluctuations is different than that of in- 
phase modulation. We found large amplitude Weiss-type 
oscillation at very low magnetic field even when band- 
width becomes non-oscillatory. For conventional 2DEG, 
effect of the out-of-phase modulation remains same as 
graphene though the bandwidth of 2DEG is oscillatory. 
This high amplitude Weiss-type oscillation is due to the 
magnetic field dependent phase factor which plays an im- 
portant role here, unlike the Weiss oscillation in electrical 
transport properties where oscillation washed out for the 
same amplitude of energy correction due to two modula- 
tions. 
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FIG. 9: Plots of the fluctuation in the Helmholtz free energy 
vs magnetic field. 



Appendix A 



VII. SUMMARY 

We have studied the effect of magnetic modulation on 
thermodynamic properties of a graphene sheet. The re- 



Here we calculate the asymptotic expression of the 
DOS of graphene in presence of a modulated magnetic 
field: B = [B + B\ cos (qx)}z. The total energy upto the 



fajO g y2n+G n cos (qxo), where 



first-order in V£ is E^'™ 



G n = l/9^7^ e -"/ 2 [L n -i( M ) - L n (u)] and u = q 2 l 2 /2. 
Using e~ u l 2 L n {u) ~ (ny/nu)~ 1 / 2 cos[2 y/nu — 7r/4] for 
higher values of n, we get the asymptotic expression of 



Ct^i as 



where 



2VI / n \ 1/4 



( — (50 sin 



x sin 



- ) sin -^ 2\/nu > 

>,/ I 4 J 



2E 
qlohujg 

qlpfihjg 
2E 



1/2 



'7'o 



Sin {^°/i:-i} 



(Al) 



Now we use impurity broadened Landau levels in limiting 
case. The self-energy can be written as^ 



a dxg 



£-(*) = r 3 W^ 



1 



E-E° n 'Z-X-(Ey 



(A2) 



where To is the broadening of the Landau levels due to 
impurities. By determining the imaginary part of the 
self-energy we can get the DOS through 



D(E) = Im 






(A3) 



By using residue theorem we get 

, 2 f a dx E - £-(£) - G n cos(qx ) 



T,-{E) 



2irTt 



x cot 



(\/2/zw 9 ) 2 
7tLq_E f a dxo 



{E - Z- {E)-G n cos (qx Q )} 2 



,, a 

7T 



(^ 5 ) 2 Jr, 



x cot 



tt£; 



r {£-2IT(£)-2G„cos(gxo)} 



(V2?iw g ) 2 
Separating T,~(E) into real and imaginary parts as 

' sin u + i sinh t> 



E-(25) = A(25) 






coshw — cosu 



(A4) 






[E-2A(E)-2G„cos(qx )], 



and w = (I ^^ 2 T(E). In the limit of large collisional 
broadening, 7rL 3> fiw 9 , after expanding the hyperbolic 
term with respect to the small quantity exp(— v) up to 
first-order one obtains 



T(E) _ irTlE 



(fiW g )2 



1 + 2 exp 



nET 



a dx r irE 



x / cos 

/o a 



(S-2G„cos( g a;o))} 



(A5) 
After first iteration, we have Y(E)/2 = nTlE/(fuv g ) 2 , 



and then putting it back into Eq. (|A5j) and using e 
E/(hujg), we get 



r(e) ^ Trrge 

2 "" {fiUg) 


[™ HO 2 } 


X COS (7T£ 2 ) 


X |l-f 


! m (e7kJ s ) 3 sin 2 fgZ e - 7) | 



(A6) 



where 



(KSiVa 



5 c^i' (A7) 



<?/() 



a 



Using Eqs. (|A3|) and (|A6j) . we get the density of states 
D(e) as given in Eqs. (fl~5l) and (ITBT) . 
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